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Abstract. Spectral properties of the x-ray pulses, generated by perfectly uniform 
cylindrical implosions of tungsten plasma with parameters typical of wire array z- 
pinches, are investigated under the simplifying assumption that the final stage of 
the kinetic-to-radiant energy conversion is not affected by the magnetic field. The 
x-ray emission is shown to be generated within a narrow (sub-micron) radiation- 
dominated stagnation shock front with a "supercritical" amplitude. The structure of 
the stagnation shock is investigated by using two independent radiation-hydrodynamics 
codes, and by constructing an approximate analytical model. The x-ray spectra are 
calculated for two values of the plasma column mass, 0.3 mg cm~^ and 6 mg cm^^, 
with a newly developed two-dimensional radiation-hydrodynamics code RALEF-2D. 
The hard component of the spectrum (with a blackbody-fit temperature of 0.5-0.6 keV 
for the 6-mg cni^^ mass) originates from a narrow peak of the electron temperature 
inside the stagnation shock. The softer main component emerges from an extended 
halo, where the primary shock radiation is reemitted by colder layers of the imploding 
plasma. Our calculated x-ray spectrum for the 6-mg cm^^ tungsten column agrees 
well with the published Sandia experimental data (Foord et al 2004 Phys. Rev. Lett. 
93, 055002). 
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1. Introduction 

Wire array z-pinches proved to be one of the most efficient and practical way to generate 
multi-terawatt pulses of quasi-thermal x-rays with duration of a few nanoseconds [1, 2]. 
This implies reach potential for many applications, in particular, as an attractive driver 
for inertial confinement fusion (ICF) [3, 4]. Prom theoretical point of view, the wire array 
z-pinch is a complex physical phenomenon: its adequate modelling requires sophisticated 
multi-dimensional magnetohydrodynamic (MHD) simulations of a complex plasma- 
metal configuration, whose dynamics at a later stage is strongly influenced by radiative 
processes [5, 6, 7]. 

The focus of this paper is on a single specific aspect of this phenomenon, namely, on 
the physics of conversion of the kinetic energy of the imploding pinch into the emerging 
pulse of x-rays. We solve an idealized one-dimensional (ID) problem where a high- 
velocity cylindrical implosion of tungsten plasma is stopped according to the laws of pure 
radiation hydrodynamics, with a possible role of MHD effects neglected. We believe that, 
despite a seemingly oversimplified statement of the problem, thus obtained results will 
prove to be useful for interpretation and subsequent modelling of the x-ray pulses and 
their spectra in wire array z-pinches optimized for maximum x-ray power. The emission 
that wc calculate should primarily refer to the main part of the observed pulses, i.e. to 
the x-ray fiux in a narrow time window around the peak of the x-ray power with the 
full width at half maximum (FWHM) roughly equal to its rise time (~ 5 ns in a typical 
18-19 MA shot on the Z machine at Sandia [8]). Note that thus defined main pulse may 
only contain about 50% of the total x-ray emitted energy [8]. 

Our analysis is based on the assumption that practically all the energy radiated 
in the main pulse originates from the kinetic energy of the imploding plasma. Such 
a premise is corroborated by the latest 3D MHD simulations of imploding wire arrays 
[7]. To simplify the treatment, wc do not consider the acceleration stage of the plasma 
implosion and start with an initial state at maximum implosion velocity. Hence, we 
do not have to consider the j x B force which drives the implosion (i.e. accelerates the 
plasma but generates little entropy) because the entire kinetic energy of the implosion 
can be simply prescribed at the initial state. The resistive (Ohmic) dissipation of the 
electromagnetic energy was found to be neghgible anyway [7]. More difficult is to 
justify our neglect of the MHD effects at the final stage of plasma deceleration and 
x-ray generation. Although some additional arguments to this point are given in the 
summary, we must admit that full clarification of this issue remains for future work. 

Under the assumptions made, we find that the kinetic energy of the implosion is 
converted into radiation when the plasma passes through a stagnation shock near the 
axis. In wire arrays with powerful x-ray emission the stagnation shock falls into the 
category of "supercritical" radiation-dominated (RD) shock fronts [9]. Such fronts are 
characterized by an extremely small width, and their thermal structure is primarily 
controlled by emission and transport of thermal radiation. Our results imply that 
adequate modelling of the temperature and density profiles across the narrow stagnation 
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shock front is key to understanding the x-ray spectra emitted by powerful wire array 
z-pinches. 

The paper is organized as follows. In section 2 we describe the initial state of 
the simulated plasma flow; section 3 elaborates on the two numerical codes used for 
simulations. Section 4 is devoted to the detailed analysis of the stagnation shock 
structure: we construct an approximate analytical model, which is corroborated by 
numerical simulations and allows simple evaluation of the key plasma parameters in the 
shock front. In section 5 we present the calculated spectra of the x-ray pulses, radial 
profiles of the spectral optical depth, spectral x-ray images of the plasma column. The 
emergent x-ray spectra have been calculated for two values of the linear pinch mass, 
0.3 mg cm~^ and 6 mg cm^^, by employing a newly developed two-dimensional (2D) 
radiation-hydrodynamics code RALEF-2D. 



2. Initial configuration 

We choose the simplest initial configuration that allows reproduction of the principal 
characteristics of the observed x-ray pulses. We start with a cyhndrical shell of 
tungsten plasma, converging onto the pinch axis r = with an initial implosion velocity 
^(0, r) = — f/o {Uq > 0) constant over the shell mass. The imploding shell is supposed to 
have sharp boundaries at r = ri{t) and r = ri(t) + Aq. Once the radial velocity peaks at 
f/o, the implosion can be treated as "cold" in the sense that the plasma internal energy 
is small compared to its kinetic energy, the role of pressure forces is negligible, and the 
shell thickness freezes at a constant value Aq. We begin our simulations at time t — 
when the inner shell edge arrives upon the axis, i.e. when ri = '"i(O) = 0. 

The density distribution across the imploding shell is assumed to have been uniform 
at earlier times, when the inner shell radius was ri{t) ^ Aq. In a cold implosion this 
leads to the initial radial density profile of the form 




where is the linear (per unit cylinder length) mass of the shell. 

In this paper we present simulations for two cases, namely, case A (referring to the 
5-MA Angara-5-1 machine in Troitsk, Russia) and case Z (referring to the 20- MA Z 
accelerator at Sandia, USA). In both cases we used the same values of the implosion 
velocity and shell thickness, 

Uo = 400 km s"^ = 4 X 10^ cm s-\ Aq = 2 mm. (2) 

The peak implosion velocity of 400 km s^^ has been inferred from the experimental data 
for optimized shots on both the Angara-5-1 [10] and the Z machines [11], and confirmed 
by numerical simulations [5, 7]. For the given Uq, the shell thickness Aq is set equal to 
2 mm to conform with the observed x-ray pulse duration of 5 ns (FWHM) [1, 2, 12, 13]. 
Note that if, in addition, we assumed a 100% instantaneous conversion of the kinetic 
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energy into x-rays, we would obtain a rectangular x-ray pulse of duration 

to = ^ = 5 ns (3) 
with the top nominal power 

Po-^ (4) 
° 2Ao ^ ^ 

(per unit cylinder length). 

Thus, the only parameter that differs between the cases A and Z is the linear mass 

mo of the imploding shell. In simulations wc used the values 

f 0.3 mg cm~\ case A, , . 

^ 1 n -1 7 ^ 

[ o.U mg cm , case /, 

which are representative of a series of optimized (with respect to the peak power and 
total energy of the x-ray pulse) experiments at a 3 MA current level on Angara-5-1 
[13], and at a 19 MA current level on Z [8]. These two values of mo correspond to the 
nominal powers 

_ I 4.8 TW cm-\ case A, 
° ~ \ 96 TWcm-i, case Z, ^ ^ 

which are close to the peak x-ray powers measured in the corresponding experiments. 

The final parameter needed to fully specify the initial state is the initial temperature 
To of the imploding plasma. In both cases we used the same value Tq — 20 eV, which 
falls in the 10-30 eV range inferred from the theory of plasma ablation in multi-wire 
arrays [14, 15, 16] and confirmed by direct MHD simulations of the wire-corona plasma 
[5]. The sound velocity in a 20-eV tungsten plasma, Cg ~ (0.5-1.0) x 10® cm s^^, implies 
implosion Mach numbers as high as Uq/cs ~ 40-80 — which fully justifies the above 
assumption of a cold implosion. 



3. The DEIRA and the RALEF-2D codes 



Numerical simulations have been performed with two numerical codes that are based 
on different numerical techniques and include fully independent models of all physical 
processes, namely, with a ID three-temperature (3T) code DEIRA [17, 18], and a 2D 
radiation- hydrodynamics code RALEF-2D [19]. Because of strongly differing physical 
models and numerical capabilities, the results obtained with these two codes are to a 
large extent complimentary to one another. The 2D RALEF code was used to simulate 
our ID problem simply because we had no adequate ID code with spectral radiation 
transport at hand. 



3.1. The DEIRA code 

The ID 3T DEIRA code was originally written to simulate ICF targets [17]. It is 
based on one-fiuid Lagrangian hydrodynamics with a tensor version of the Richtmyer 
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artificial viscosity and different electron, Tg, and ion, 7], temperatures. Included also are 
the electron and the ion thermal conduction, as well as the ion physical viscosity. The 
model for the electron and the ion conduction coefficients is based on the Spitzer formula, 
modified in such a way as to match the experimental data near normal conditions [18]. 
The equation of state is based on the average ion model [20], which accounts for both 
the thermal and the pressure ionization at high temperatures and/or densities, as well 
as for realistic properties of materials near normal conditions. 

Energy transport by thermal radiation is described by a separate diffusion equation 
for the radiative energy density per = asT^, expressed in terms of a separate radiation 
temperature Tj-; here as is the Stefan constant. The energy relaxation between the 
electrons and the radiation is expressed in terms of the Planckian mean absorption 
coefficient kp\, while the radiation diffusion coefficient is coupled to the Rosseland 
mean /cr. The absorption coefficients kp\ and /cr are evaluated in-line as a combined 
contribution from the free-free, bound-free and bound-bound electron transitions 
plus the Thomson scattering. For the bound-bound and bound-free transitions, an 
approximate model, based on the sum rule for the dipole oscillator strengths [21], is 
used. The subset of the three energy equations is solved in a fully implicit manner by 
linearizing with respect to the three unknown temperatures T^, 7], and T^.. 

3.2. The RALEF-2D code 

RALEF-2D (Radiative Arbitrary Lagrangian-Eulerian Fluid dynamics in two 
Dimensions) is a new radiation-hydrodynamics code, whose development is still 

underway [19]. Its hydrodynamics module is based on the upgraded version of the 
CAVEAT hydrodynamics package [22]. The one- fluid one-temperature hydrodynamic 
equations are solved in two spatial dimensions [in either Cartesian (x, y) or axisymmetric 
(r, z) coordinates] on a multi-block structured quadrilateral grid by a second-order 
Godunov-type numerical method. An important ingredient is the rezoning-remapping 
algorithm within the Arbitrary Lagrangian-Eulerian (ALE) approach to numerical 
hydrodynamics. The original mesh rezoning scheme, based on the Winslow equipotential 
method [23], proved to be quite efficient for the interior of the computational domain, 
if the mesh is smooth along the boundaries; in RALEF, a new high-order method for 
rezoning block boundaries has been implemented to this end. 

New numerical algorithms for thermal conduction and radiation transport have 
been developed within the unified symmetric semi-implicit approach [24] with respect 
to time discretization. The algorithm for thermal conduction is a conservative, second- 
order accurate symmetric scheme on a 9-point stencil [25]. Radiation energy transport 
is described by the quasi-static transfer equation 



for the spectral radiation intensity Jj^ = /,^(t,x, il); the term c ^dl^/dt, where c is the 
speed of light, is neglected. A non-trivial issue for spatial discretization of equation (7) 




(7) 
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together with the radiative heating term 

oo 

Q^ = - div J du J I^rtdn (8) 

47r 

in the hydrodynamic energy equation, is correct reproduction of the diffusion hmit on 
distorted non-orthogonal grids [26]. In our scheme, we use the classical Sn method to 
treat the angular dependence of the radiation intensity Iu(t,yi,fl), and the method 
of short characteristics [27] to integrate equation (7). The latter has a decisive 
advantage that every grid cell automatically receives the same number of light rays. 
Correct transition to the diffusion limit is achieved by special combination of the 
first- and second-order interpolation schemes in the finite-difference approximations to 
equations (7) and (8). More details on the numerical scheme for radiation transfer are 
to be pubhshed elsewhere. 

In the present work we used the equation of state, thermal conductivity and spectral 
opacities provided by the THERMOS code [28], which has been developed at the Keldysh 
Institute of Applied Mathematics (Moscow). The spectral opacities are generated 
by solving the Hatree-Fock-Slater equations for plasma ions under the assumption of 
equilibrium level population. In combination with the equilibrium Planckian intensity 
B^, used in (7) as the source function, the latter means that we treat radiation transport 
in the approximation of local thermodynamic equilibrium (LTE) — which is justified 
for relatively dense and optically thick plasmas considered here. 

] T=0.25 keV, p=0.01 g/cc \ 




, , — ,— ■-, i i — I I I I 

0.1 1 10 

Photon energy hv (keV) 



Figure 1. (Colour online) Spectral absorption coefBcient k^, [cm^^] of tungsten at 
p = 0.01 g cni^^, T = 250 eV used in the present simulations: shown are the original 
data from the THERMOS code (thin solid curve) together with the group-averaged 
values for 8 (dashed) and 32 (thick solid) selected spectral groups. 

The transfer equation (7) is solved numerically for a selected number of discrete 
spectral groups [i/j,i/j+i], with the original THERMOS absorption coefficients /c^ 
averaged inside each group j by using the Planckian weight function. Two different 
sets of frequency groups are prepared for each code run: the primary set with a smaller 
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number of groups (either 8 or 32 in the present simulations) is used at every time step in 
a joint loop with the hydrodynamic module, while the secondary (diagnostics) set with 
a larger number of groups (200 in the present simulations) is used in the post-processor 
regime at selected times to generate the desired spectral output data. An example of 
the spectral dependence of /c^, provided by the THERMOS code for a tungsten plasma 
at p = 0.01 g cm~^, T — 250 eV, is shown in figure 1 together with the corresponding 
group averages used in the RALEF simulations. 

3.3. Numerical setup for the RALEF simulations 

To test the sensitivity of the results with respect to spectral radiation transport, we did 
our simulations with two selections for the primary set of frequency groups, namely, 

• with 8 groups delimited by the photon energies 

huj = 10"^ 0.1, 0.2, 0.4, 0.8, 1.5, 3.0, 6.0, 10.0 keV, (9) 

• and with 32 groups delimited by 

huj = 10-^ 0.02, 0.04, 0.07, 0.1, 0.119, 0.1414, 

0.168, 0.20, 0.238, 0.2828, 0.336, 0.4, 0.476, 

0.5656, 0.672, 0.8, 0.952, 1.1312, 1.344, 1.5, 

1.785, 2.121, 2.52, 3.0, 3.57, 4.242, 5.04, 

6.0, 7.14, 8.484, 10.08, 12.0 keV. (10) 

The 200 spectral groups of the secondary (diagnostics) frequency set were equally spaced 
along \n(hu) between hui — 0.01 keV and hu2oi — 10 keV. The angular dependence of 
the radiation intensity was calculated with the 5'i4 method, which offers 28 discrete ray 
directions per octant. 

The simulated region occupied one quadrant < (f) < 90° of the azimuth angle 
(f) with reflective boundaries along the x- and y-axes. Near the geometrical centre 
X = y = 0, a rigid transparent wall was placed at r = tq = 10 fim with the boundary 
conditions of u{t, r^) = and zero thermal flux. Thermal radiation passed freely through 
this cylindrical wall and was reflected by the two perpendicular reflective boundaries. 
Two variants of the initial polar mesh were used: a x nr = 50 x 250 mesh in case A, and 
a X rij. = 60 X 600 mesh in case Z. At the outer boundary (initially at r — Rq — 2 mm), 
the boundary conditions of zero external pressure and zero incident radiation flux were 
applied. 

4. Stagnation shock 

4.1. General picture 

Upon arrival at the axis, the imploding plasma comes to a halt passing through a 
stagnation shock. In our situation the speciflc nature of this shock is defined by the 
dominant role of the radiant energy exchange. A detailed general analysis of the 
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structure of such RD shock waves is given in [9, Ch. VII]. Perhaps the most salient 
feature of an RD shock with a supercritical amplitude is a very narrow local peak of 
matter temperature immediately behind the density jump [9]. This temperature peak 
manifests itself as a hair-thin bright circle at r = rg ~ 43 fim on the 2D temperature 
plot in figure 2. 




X (mm) X (mm) 

Figure 2. (Colour) 2D contour maps of matter temperature at t = 3 ns in case A 
(as calculated with the RALEF code): frame (b) is a blow up of the central part of 
the full view (a). Colour represents matter temperature T. A thin dark- red circle at 
r — 4:3 /im marks the position of the stagnation shock. 

To achieve an adequate numerical resolution of the RD shock front, one needs a 
very fine grid that can only be afforded in ID simulations. Figures 3 and 4 show the 
density and temperature profiles across the stagnation shock at t = 3 ns as calculated 
with the ID DEIRA code on a uniform Lagrangian mesh with 20 000 mass intervals. If 
we define the shock-front width Arg to be the FWHM of the hump on the Tg profile, 
we obtain Avg = 0.5 fim in case A, and Arg = 0.3 /im in case Z. The peak values of the 
electron temperature are calculated to be Tgp = 0.35 keV in case A, and Tgp = 0.54 keV 
in case Z. 

In a model where one distinguishes between the electron and the ion temperatures 
but ignores viscosity and ion heat conduction, the shock front has a discontinuity in the 
density and Ti profiles. In figures 3 and 4 this discontinuity is smeared out (roughly over 
3 mesh cells) by the artificial Richtmyer-type viscosity, present in the DEIRA code. The 
electron temperature Tg, which exhibits a prominent hump over the virtually constant 
radiation temperature T^, is continuous because the electron thermal conduction plays 
a significant role. Clearly, the plasma inside the dense part of the Tg hump must be 
intensely loosing energy via thermal radiation. 

Figures 3 and 4 demonstrate prominent peaks of the ion temperature Ti, whose 
maximum values Tjp ^ 20 keV significantly exceed the peak electron temperature Tgp. 
This fact, however, turns out to be rather insignificant for the radiative properties of the 
imploding plasma. Indeed, if we assume that the kinetic energy of the infalling plasma is 
fully converted into the ion thermal energy within a density-temperature discontinuity 
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Figure 3. (Colour online) Density and temperature profiles across the stagnation 
shock in case A as calculated with the DEIRA code for t = 3 ns. The effective width 
Ars of the shock front is defined as the FWHM of the local peak of the electron 
temperature To. The velocity profile in the shock frame can be easily restored from 
the density plot and the condition pv — constant, which is quite accurately observed 
in the displayed region. 
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Figure 4. (Colour online) Same as figure 3 but for case Z. 



and ignore the plasma preheating before the shock, we calculate an after-shock ion 
temperature of 

Ti+ = = = 100 keV; (11) 

here 7 = 5/3 is the adiabatic index of the ideal gas of plasma ions, m; is the mass of a 
tungsten ion. The DEIRA simulations demonstrate much lower peak ion temperatures 
because the preheating of the pre-shock plasma electrons, followed by their adiabatic 
compression in the density jump, consumes a large portion of the initial ion kinetic 
energy (in a collisionless manner via ambipolar electric fields). As a consequence, even 
before the coUisional electron-ion relaxation sets in, the post-shock electrons with a 
temperature of ^ 0.4 keV already contain almost twice as much energy as the 
post-shock ions with a temperature of T; ^ 20 keV. Hence, the subsequent coUisional 
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electron-ion relaxation does not significantly affect the Tg profile. This fact has also 
been verified directly: having performed additional DEIRA runs in the 2T mode (i.e. 
assuming = T\ = T), we obtained T profiles that were hardly distinguishable from 
the Tc profiles in figures 3 and 4 (the difference between the peak values Tp and Tep did 
not exceed 3%). Thus, the approximation of a single matter temperature T = = Ti, 
used in the 2D RALEF code, is well justified for our problem. 



4-2. Analytical model 

The theory of RD shock fronts, developed by Yu. P. Raizer [29] and described in his 
book with Ya. B. Zel'dovich [9], applies to planar shock waves in an infinite medium, 
which eventually absorbs all the emitted photons. We, in contrast, are dealing with 
a finite plasma mass, which lets out practically all the radiation fiux generated at the 
shock front. In addition, the electron heat conduction, ignored in Raizer's treatment, 
plays an important role in formation of the temperature profile across the shock front. 
Hence, we have to reconsider certain key aspects of the Raizer's theory in order to obtain 
an adequate model for the stagnation shock in imploding z-pinch plasma. 



4-2.1. General relationships To construct an analytical model of the plasma fiow, 
we have to make certain simplifying assumptions. First of all, we assume a single 
temperature T for ions and electrons and employ the ideal-gas equation of state in the 
form 

p = ApT, e=^T, (12) 
7-1 

where p is the pressure, e is the mass-specific internal energy, and A and 7 > 1 are 
constants. The thermodynamic properties of the tungsten plasma in the relevant range 
of temperatures and densities arc reasonably well reproduced with 

' 13 MJ g-^ keV"^ in case A, 
20 MJ g-^ keV"^ in case Z, 



. I ±0 ivid g, V m uase i\, (■\n\ 

^ ^ 1 on ^^T ^-1 ;„ „„„„ rj (l>^j 



{1.29 in case A, ,^ 
1.33 m case /. 

At each time t the entire plasma fiow can be divided into three zones: the inner 
stagnation zone (the compressed core) at < r < rg behind the shock front, the 
shock front itself confined to a narrow layer around r = rg, and the outer layer of 
the unshocked inf ailing material at r > rg. In the stagnation zone the plasma velocity 
is small compared to C/q, while the temperature and density are practically uniform and 
have the final post-shock values olT — Ti — Ti{t), p — pi — pi{t). Because the plasma 
fiow in the stagnation zone is subsonic, pressure is rapidly equalized by hydrodynamics, 
while temperature is equalized by efficient radiative heat conduction; note that typical 
values of the mean Rosseland optical thickness of the stagnant core lie in the range 
Tc^Ros — 10-100. Numerical simulations confirm that spatial density and temperature 
variations across the compressed core do not exceed a few percent. 
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We identify the shock radius rg = rs{t) with the density (and velocity) discontinuity, 
which is always present in sufficiently strong shocks once a zero physical viscosity is 
assumed [9]. Here and below the term "shock front" is applied to a narrow layer with 
an effective width of Arg ^ rg, where the matter (electron) temperature T exhibits a 
noticeable hump above the radiation temperature T,-; see figures 3 and 4. Outside the 
shock front at r > Tg lies a broad preheating zone, which extends virtually over the 
entire unshocked material and has a width well in excess of the shock radius Tg. In this 
region the infalling plasma is preheated due to interaction with the outgoing radiation; 
in the process, it is also partially decelerated and compressed. 

To analyze the structure of the shock front, we use the three basic conservation 

laws 

pv = —j — constant, (15) 
p + pv^ — constant, (16) 

= constant, (17) 

governing a steady-state hydrodynamic flow without viscosity across a planar shock 
front [9]; here 

^ = e + ^ = ^T (18) 
P 7-1 

is the specific enthalpy, Sc and Sj- arc the energy fluxes due, respectively, to the electron 
thermal conduction and radiation transport. Equations (15)-(17) are written in the 
reference frame comoving with the shock front: in this frame the plasma velocity f < 0. 
To avoid confusion, we use symbol "f " for the plasma velocity in the shock frame, and 
symbol "li" for the plasma velocity in the laboratory frame. The velocity of the shock 
front in the laboratory frame is Ug — drg/di. Clearly, the planar conservation laws (15)- 
(17) can be applied over a narrow front zone with |r — rg] <^ Tg but not across the broad 
preheating zone, where the effects of cylindrical convergence become significant. 




4-2.2. Parameters of the stagnant core We begin by deriving a system of equations, 
from which the parameters of the stagnant plasma core can be evaluated. Although 
the sought-for quantities formally depend on time t, time appears only as a parameter 
in the final equations. The final post-shock plasma state can be determined in the 
approximation of zero thermal conduction, which redistributes energy only locally, in 
the immediate vicinity of the shock front. Without thermal conduction, the density and 
the temperature should have profiles shown qualitatively in figure 5: the density jump 
from p = p_ to p = p+ is accompanied by the jump in temperature from T — to 
T — T^. Immediately behind (in the downstream direction) the temperature peak r+ 
lies a narrow relaxation zone to the final state (pi, Ti), where the excess thermal energy 
between the T+ and Ti states is rapidly radiated away. 

As was rigorously proven by Ya. B. Zel'dovich [30], the preheating temperature 
TL at the entrance into the density jump can never exceed Ti. How does T_ compare 
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Figure 5. (Colour online) Schematic view of the density, p, and temperature, T, 
profiles across an RD shock front in the approximation of zero viscosity and heat 
conduction. The relaxation zones before and after the density jump are due to energy 
transport by radiation. 



with Ti, depends on whether the RD shock is subcritical or supercriticaL A critical 
amphtude of an RD shock front corresponds to the condition [29] 

ATi 

aT^ = piMsCi ^ poUo (19) 

7 - 1 

that the one-sided radiation energy flux aT^ becomes comparable to the hydrodynamic 
energy flux piUgCi behind the shock front; here a is the Stefan-Boltzmann constant. In 
our case the RD shock becomes supercritical when the post-shock temperature exceeds 
the critical value of 

ri,e.^0.3p;/^/^^keV, (20) 

where Po,g/cc is the initial density (1) at r = rg in g cm^^. A posteriori, having 
calculated Ti and rg from the equations given below, we verify that our shock fronts 
are supercritical. As shown by Yu. P. Raizer [29], when a supercritical RD shock wave 
propagates in an infinite medium, it has T_ ^ Ti. If the optical thickness of the 
unshocked material is not exceedingly small, this is also true in the case of a finite 
plasma size; the latter applies to all configurations considered in this paper and is 
directly confirmed by the profiles in figures 3 and 4. 

To avoid treatment of the non-planar preheating zone, we make a simplifying 
assumption that partial deceleration of the infalling plasma in the preheating zone can 
be neglected, i.e. that one can set p_ = po and f_ = —{us + Uq), where po is calculated 
from (1) at r = r^, and f _ is the plasma velocity at the entrance into the jump in the 
shock front frame. As is demonstrated in §16 of chapter VII in [9], for 7 — 1^1 this 
approximation is accurate to the second order with respect to the small parameter 

r/ioo = (7-l)/(7 + l)- (21) 
In fact, it has already been used in condition (19). 
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Now we can apply equations (15), (16) of mass and momentum balance between 
the states p_,TL and pi,Ti: 

Po{us + C/o) = Pikil = j, (22) 
p- + Po{us + Uof ^ pi + pivf; (23) 

here po — Po('^s); and Vi is the unknown material velocity behind the shock front. 
Because the shock wave propagates over a falling density profile po{r) oc r~^, a uniform 
density distribution behind the front implies that the post-shock density pi{t) decreases 
with time, and material in the stagnation zone expands. The expansion velocity u is 
small compared with Uq but not with the front velocity Us- As a consequence, we cannot 
simply put Vi = —Us. 

By virtue of (22) and (12), equation (23) can be transformed to 

Vi{l-Vi){u. + Uor = ^(l-Vi^), (24) 

Pi V TiJ 

where 

m = — = ——jr 25) 

Pi Us + Uq 

is the inverse of the compression factor. Restricting our treatment to the case of 
supercritical RD shocks, where T_ 71, we get 

/IT, 

^ <^^' 

In our model r]i is a small parameter, which is even smaller than the inverse compression 
factor r^ioo in the infinite-media. Keeping this in mind, in all the algebra below we 
consistently retain only the zeroth and the first terms with respect to this parameter. 

A subtle point here is that we cannot directly use equation (17) of the energy balance 
across the shock front. Quasi- uniform density and temperature profiles in the stagnation 
core ensue from the rapid redistribution of thermal energy over the entire mass of this 
zone by means of radiation. Hence, the post-shock thermal energy calculated from the 
local condition (17) may differ considerably from the required average value. To obtain 
the latter, we use the condition of global energy balance. For a similar reason, we employ 
the equation of global mass balance to establish the relationship between the radius Tg 
and the velocity Ms of the shock front. 

The total mass nis of the compressed core can be expressed as 

mg = ms(i) = Trr^pi = (r-s + t/o^) , (27) 

^0 

which, by virtue of (25) and (1), yields 



2rj 



^ Uot. (28) 



1 - 2?7i 

Since rji varies only slowly with time (this can be verified a posteriori), equation (28) 
implies 

_ drs 2r/i 

Us = — = — Uq, Ts = Ust. (29) 

at 1 — zfji 
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Combining equations (29) and (26) and omitting the second and higher order terms 
with respect to rji, we obtain 

ATi Us _ 2ATi 

f/^ + 4ATi' TFo ~ Ui + 2ATi ■ 
The global energy balance for the imploding plasma mass can be expressed as 
d 



= ^ = ^^^^ /orvN 



nisei + (mo - "^s)^ 



= 0, (31) 



where Px is the total (per unit cylinder length) power of x-ray emission, which escapes 
through the outer boundary. If Px = 0, we obtain a simple "conservative" result 

ei^lu^, (32) 



which yields 



rr. _rj. _ 1''^ j r2 _ / 1-8 kcV lu casc A, 

- - ^ - 1.3 keV in case Z. ^^^^ 



When equation (32) is used with the realistic equation of state of tungsten, provided 
by the THERMOS code, it yields T^^ ^ 0.8 keV in case A, and T^^ ^ 0.95 keV 
in case Z. It is this post-shock temperature that one would calculate, having literally 
applied the Raizer's model to a planar stagnation shock in an infinite medium. In 

our non-conservative situation, where most of the radiation flux escapes the imploding 
plasma, the final post-shock temperature Ti is significantly lower than Tioo- 

To close up our analytical model, we need an expression for Px- If we assume the 
unshocked infalling plasma to be transparent for the outgoing radiation, we can write 

Px^2nrsaT^, (34) 

which means that the opaque compressed core of radius Vs radiates as a black body with 
a surface temperature Ti. Clearly, such a situation should correspond to sufficiently 
small values of mg, and our case A, as will be seen below, falls into this category. 

An additional approximation that we make when opening the brackets in (31) is 
neglect of the term mgdei/dt compared to eidrng/dt: this spares us the need to solve a 
differential equation with practically no loss of accuracy. As a result, upon substitution 
of (27), (29), (30 and (34) into (31), we arrive at the following equation for determination 
oiTi^Tu^Tuit) 

KmoJ U^ + 4ATu 

where Tioo is given by (33). Here we introduced a separate notation Ti* for the post- 
shock temperature Ti, calculated from (35) in the optically thin approximation for the 
pre-shock plasma, when expression (34) is applied. Having found Ti = Tu from (35), 
we calculate rji, Us, Vs and pi from equations (30), (29) and (25), respectively, and this 
completes our analytical model for the plasma parameters in the stagnation core. 

Figure 6 compares the values of Ti = Ti*, calculated from equation (35), with those 
obtained in the DEIRA and RALEF simulations. A very good agreement is observed in 



Tioo, (35) 
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Figure 6. (Colour online) Time dependence of the post-shock temperature Ti: 
solution of equation (35) Ti = Ti* — Ti*(t) (solid curves) is compared with the results 
of the DEIRA (open circles) and the RALEF (crosses) simulations. 



case A, where the optical thickness of the pre-shock plasma at different frequencies has 
moderate values around 1 (see figure 16 below). The agreement becomes worse in case Z, 
where Tg reaches values around 10 and higher (see figure 19 below): then equation (34) 
significantly overestimates the radiative energy loss. From the above analysis it follows 
that the true value of Ti should be in the range Ti* < Ti < Tioo; when Ts > 1 increases, 
the difference Ti — Tu grows and Ti approaches the limiting value of Ti = Tioo- Note 
that, when considered as a function of the total imploding mass mo at a fixed value 
of Uq, the post-shock temperature Ti grows with mo firstly because Ti^ increases [as it 
follows from (35)], and, secondly, because the difference Ti — Ti* > becomes larger for 

Ts > 1. 

4-2.3. Temperature peak in the shock front In addition to the post-shock parameters, 
one would like to have an estimate for the peak matter (electron) temperature Tp inside 
the shock front (see figures 3 and 4), which defines the hard component of the emitted x- 
ray spectrum. Such an estimate, however, cannot be obtained without a proper account 
for thermal conduction. With the conduction energy flux given by 



where k is the conduction coefficient, the temperature T becomes a continuous function 
across the density jump, while Se is discontinuous; the radiation energy flux S^, on the 
contrary, is everywhere continuous [9]. A qualitative view of the density, temperature 
and the energy flux profiles across a supercritical RD shock front with strong thermal 
conduction is shown in figure 7. 

In the shock structure shown in figure 7 one can identify four hydrodynamic states: 
state at the foot of the conduction-preheated layer before the density jump, state "— " 
at the entrance into the density jump, state upon the exit from the density jump, 
and state 1 behind the post-shock relaxation zone. The effective width (FWHM) of the 
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Figure 7. (Colour online) Schematic view of the structure of a supercritical RD shock 
front with non-zero thermal conduction and zero viscosity. Shown are the density, p, 
and the temperature, T, profiles (a), as well as the profiles of the conductive, So, and 
the radiative, ^r, energy fiuxes (b). At r — Vg the density p and the conductive fiux 
are discontinuous. 



conduction-preheated layer is h_; the effective width of the post-shock relaxation layer 
is the effective width of the entire shock front is the sum of the two, 

Ars = h^ + h+. (37) 

In our case h- is determined primarily by thermal conduction, whereas for /i+ both 
the radiant emissivity and thermal conduction are important. Because /i_ is much 
shorter than the shock radius Ts, we can assume that the state lies at the end of 
the broad radiation-preheating zone and, as in the previous subsection, ignore partial 
plasma compression and deceleration in the latter. Then, the plasma parameters in the 
four mentioned states can be represented as in table 1. The parameters in states and 
1 are known from the previous subsection. Here we have to evaluate /i+ and Tp. 

We derive a system of approximate equations for the three unknowns h+ and 
Tp by successively applying the general equations (16), (17) of momentum and energy 
balance three times, for transitions between states and "— ", between states "— " and 
and between states "+" and 1. For each of the three transitions we define a 
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Table 1. Plasma parameters at four characteristic states inside the RD shock front. 
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corresponding inverse compression factor 

_ Po _ P- _ P+ /OQN 

77_ = — , = — , 77+ = — ; (38) 

P- P+ Pi 

evidently, we must have 

77_77s77+ = 771. (39) 

For the first transition between states and "— " we can neglect the coupling 
between radiation and matter because here the plasma density p ~ po is low compared 
to the compressed state. The latter means that S-cq ~ 5'r_, and equations (16), (17) can 
be written as 



A(rp-Tir7_)=77_(l-77_)^o', (40) 

% = -^(7^p-7^i)-|(l-r/!). (41) 

J 7 - 1 2 

Because Tp — Ti < Ti, we have 1 — 77. 1, which enables us to reduce equations (40) 
and (41) to 

l-rj^^A{T^-T,)/vl (42) 

^-;^(Tp-TO. (43) 

The second transition is an isothermal density jump between states "— " and "+", 
where S^- is continuous and jumps from S^- to 8^+- Here equations (16), (17) take 
the form 

ATp = ri.vl = 7]sV-vl (44) 

% = %-%l-^^ (45) 
J J ^ 

Neglecting the second and higher order terms with respect to the small parameters 77s 
and 1 — 77- in (45), we find 

%«-^(r,-r.)-|. (46) 

J 7 - 1 2 
The third transition from state "+" to state 1 occurs in the compressed state, where 
the plasma emissivity (roughly proportional to the density p) is high, and we have to 
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account for variation of the radiation energy flux Sj-. Hence, equations (16), (17) take 
the form 

A(Ti - ij+Tp) = r,X{rt+ - 1), (47) 
= 73i(^'p - r.) + 5 ,14 {'U' - 1) - (48) 

Retaining just the leading terms with respect to the small parameter r)i, we obtain 

r)+^T,/T^ ^ piT,^p+Tp, (49) 

Sr+-S,,^^jv',^^jUl (50) 

As a final step, we express the heat conduction fluxes in terms of the corresponding 
temperature gradients, 

and the radiation flux increment 

- 5rl ~ ^ <7kpih+ (t^ - Tt) (52) 

in terms of the post-shock plasma emissivity; in (51) and arc, respectively, 
the conduction coefficients in states "— " and "+"; in (52) k-p\ is the Planckian mean 
absorption coefficient of radiation in state "+". Expression (52) is an approximation 
to the emission power of an optically thin planar layer, which is valid in both limits of 
Tp ^ Tj. = Ti and Tp ^ Tj. = Ti; the factor |(T instead of 4cr takes into account that 
h+ is the halfwidth of the T rather than the profile. From (43), (46), (50)-(52) we 
obtain the following system of three equations for evaluation of /i+ and T^: 

/i_ = (7-l)fi;_/(2jA), (53) 

1 JUI 

16akpiiT^-Tt) 



jUi ^^^2^_ ^^^^ 



Tp - Ti h+ 7-1 
For numerical estimates we use power-law approximations 

0.15Tp2keV TW cm-1 keV'S (56) 
A;pi ~ 700 ^ 700 P^^^l^f^'^^^-^ cm-\ (57) 

-tp.keV -'p.keV 

to the THERMOS data for tungsten in the relevant parameter range. 

In table 2 the analytically evaluated shock parameters are compared with those 
obtained in the RALEF simulations for t = 3 ns. Generally, the analytical model tends 
to produce higher values of the peak temperature Tp than the DEIRA and the RALEF 
codes because of an assumed sharp angle in the temperature profile (see figure 7), which 
is smeared either by artificial viscosity in the DEIRA code, or by insufficient spatial 
resolution in the RALEF simulations. It is clearly seen that, as one passes from case A 
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Table 2. Comparison of the analytically evaluated stagnation shock parameters with 
the RALEF results. 
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Tp (keV) 
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h- (/xm) 


1.08 




0.24 




h+ (urn) 


0.14 




0.08 





to a more powerful case Z, the stagnation shock becomes significantly hotter and more 
narrow — in full agreement with general properties of supercritical RD shocks [9]. 
Because of the intricate coupling between thermal conduction and radiation emission, 
no universal power-law scaling for Tp and Arg can be deduced from (53)-(57). 

5. X-ray pulse 

The 3T model of the DEIRA code is reasonably adequate for calculating the total power 
profile of the x-ray pulse (see figures 8 and 9 below), but can provide no information on 
its spectral characteristics. For this one has to solve the equation of spectral radiation 
transfer together with the hydrodynamics equations, and that is where we employ the 
RALEF-2D code. 

5.1. Power profile 

Figures 8 and 9 show the temporal x-ray power profiles Px = Px{t) as calculated with 
the DEIRA and the RALEF codes, which agree fairly well with one another, especially 
in case A. These profiles demonstrate a clear quasi-steady phase, which lasts about 4 ns 
in case A, and about 2.5 ns in case Z; at this phase Px is close to the nominal power Pq- 
A marked difference between cases A and Z is a later (by ~ 1 ns) rise of the x-ray power 
in case Z. This delay occurs because in case Z the radiation heat wave has to propagate 
through a more massive and optically thick layer of cold plasma before it breaks out 
to the surface. The overall efficiency of conversion of the initial kinetic energy into 
radiation (by t = 6 ns) is 92% in case A and 78% in case Z according to the RALEF 
data, and 94% in case A and 81% in case Z according to the DEIRA results. 

5.2. Shock structure in the RALEF simulations 

Radial density and temperature profiles in the imploding plasma, obtained with the 
RALEF code, are shown in figures 10 and 11. Despite quite different physical models, 
the RALEF and the DEIRA results agree almost perfectly in case A: we calculate 
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Figure 8. (Colour online) Temporal profile of the total x-ray emission power Px 
in case A: the results of three different numerical simulations are compared among 
themselves and with the nominal power profile, which corresponds to an instantaneous 
100% conversion of the plasma kinetic energy into x-ray emission. 
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Figure 9. (Colour online) Same as figure 8 but in case Z. 



practically the same values of the post-shock, Ti, and the peak, Tp, matter temperatures. 
Figure 10 also demonstrates that in case A the temperature peak is fairly well resolved 
in 2D simulations, although appears somewhat broader than in the ID DEIRA picture. 
Larger 2D values of the shock radius are explained by different position of the inner 
boundary (at r = 10 fim in the 2D case versus r = in the ID case). 

In case Z, on the contrary, the temperature peak is rather poorly resolved in 2D 
simulations, as one sees in figure 11 — despite a larger total number of radial mesh 
zones (600 in case Z versus 250 in case A). The reason is twofold: on the one hand, 
the temperature peak in case Z is about a factor 2 more narrow than in case A; on the 
other, a considerably larger shock radius rg causes the RALEF mesh rezoning algorithm 
to force a coarser grid along the radial direction. Nevertheless, the agreement between 
the RALEF and the DEIRA results for the post-shock, Ti, and the peak, Tp, matter 
temperatures is also fairly good. 

Radial profiles of the implosion velocity —u{r) at t = 3 ns are displayed in figure 12 
for both cases A and Z. One notices that the fluid velocity changes sign across the shock 
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Figure 10. (Colour online) Radial density and temperature profiles at t = 3 ns in 
case A obtained in the 2D RALEF simulation with 32 spectral groups. 
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Figure 11. (Colour online) Same as figure 10 but in case Z. 



front. As was already mentioned in section 4.2.2, the post-shock plasma on average 
slowly expands (i.e. has a negative implosion velocity) because the stagnation shock 
propagates over a falling density profile. Near the outer edge, the infalling plasma 
has already been significantly decelerated, especially in the more massive case Z. The 
decelerating pressure gradient is created by re-deposition of radiant energy transported 
from the stagnation shock front. 

Figure 13 shows spatial profiles of the mean ion charge Zio^ at t = 3 ns. One 
sees that tungsten ions with charges of Zion = 40-45 are present inside the stagnation 
shock front. It should be reminded here that these ion charges have been calculated 
in the LTE limit. A direct evidence that the LTE approximation is quite adequate in 
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Figure 12. (Colour online) Radial profiles of the plasma implosion velocity (minus 
the radial velocity m) at t = 3 ns in cases A and Z as calculated by the RALEF code 
with 32 spectral groups. Logarithmic scale for the radius allows to show the detailed 
structure of the compact shocked region together with the overall large scale behaviour. 
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Figure 13. (Colour online) Radial profiles of the tungsten ionization deg 

< = 3 ns in cases A and Z as calculated by the RALEF code with 32 spectral groups 

in the radiation transport module. 



our situation is a close agreement between the radiation and matter temperatures in 
figures 10 and 11. The apphcabihty of LTE can only be questioned inside the narrow 
shock front. However, the plasma density there is already so high [uc > 6 x 10^^ in 
case A, and ne > 5 x 10^^ in case Z) that non-LTE corrections to the values of z^on and 
T inside the shock front are not expected to be significant (poor spatial resolution of 
this region may, in fact, be a no less important issue). 
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5.3. Spatially integrated spectra 

5.3.1. Case A The overall x-ray spectrum emitted by the imploding pinch in case A 
at t = 3 ns is shown in two different representations in figures 14 and 15. This spectrum 
would have been observed through an imaginary slit perpendicular to the pinch axis by a 
detector without spatial resolution. More precisely, figures 14 and 15 display the spectral 
power Fy [TW cm~^ sr~^ keV~^] per unit cylinder length, obtained by integrating along 
the slit the intensity of the outgoing radiation, which propagates in direction 

O perpendicular to the pinch axis. The shown spectrum was obtained by solving the 
transfer equation (7) in the post-processor mode for 200 spectral groups of the secondary 
frequency set. In case A it turns out to be rather insensitive to the number of spectral 
groups (either 8 or 32) coupled to the hydrodynamics energy equation. 
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Figure 14. (Colour online) Spectral power of x-ray emission per unit cylinder length 
at i = 3 ns in case A: the soft part of the x-ray spectrum. The Planckian-fit curve is 
normalized to the emission peak. 
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Figure 15. (Colour online) Spectral power of x-ray emission per unit cylinder length 
t = 3 ns in case A: the hard part of the x-ray spectrum. The Planckian-fit curve is 
normalized to the hv > S keV tail of the emission. 



The plots in figures 14 and 15 demonstrate that the emitted spectrum can roughly 
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be approximated as a superposition of two Planckian curves: one with a temperature 
Tri ~ 0.11 keV, and the other with a temperature Tr2 ~ 0.34 keV. The interpretation 
of the hard component is straightforward: it is the thermal emission of the temperature 
peak Tp = T^p ~ Tr2 inside the stagnation shock. In our case this component carries 
about 16% of the total x-ray flux and is emitted by an optically thin plasma layer rather 
than by a surface of a black body. 

The soft component originates from a much broader halo around the shock front, 
at an effective radius of rem ~ 0.4 mm ^ rs = 0.043 mm. This halo is the result of 
reprocession of the original shock emission by the cold layers of the unshocked material. 
Note that the temperature Tri of the soft component is significantly lower than the 
post-shock matter temperature Ti = 0.20 keV, which implies that even in the low-mass 
case A the infalling unshocked plasma is not truly optically thin. 

Figure 16 provides more detailed information on the radial profiles of the spectral 
optical depth. It is seen that, depending on the photon energy, the optical thickness 
of the unshocked plasma can be either significantly below or significantly above unity. 
The latter means that the effective emitting layer is, in fact, not well defined, and the 
observed spectrum may exhibit significant deviations from the Planckian shape. Indeed, 
a number of prominent dips and spikes in the calculated spectrum in figures 14 and 15 
arise as a combined effect of a complex spectral dependence of the tungsten opacity, 
shown in figure 1, superimposed on a nontrivial temperature distribution inside and 
above the stagnation shock. 




Radius (mm) 

Figure 16. (Colour online) Profiles of radial optical thickness at different photon 
energies for i = 3 ns in case A. 

5.3.2. Case Z Figures 17 and 18 display the same information as figures 14 and 
15 but for a 20 times larger (6 mg cm~^) imploding mass of case Z. Here both the 
main component of the spectrum in figure 17 and the hard component in figure 18 
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correspond to roughly two times higher Planckian-fit temperatures of Tj-i = 0.21 keV 
and T^2 = 0.53 keV; the hard component carries about 7% of the total x-ray flux. 
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Figure 17. (Colour online) Same as figure 14 but for case Z. 
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Figure 18. (Colour online) Same as figure 15 but for case Z. The Planckian-fit curve 
is normalized to the hiy >6 keV tail of the emission. 

In contrast to case A, now the shock front lies at an optical depth Tg well in excess 
of unity at all frequencies: as can be seen in figure 19, at t = 3 ns the optical depth 
Ts varies in the range Xg ~ 4-100. As a result, the calculated spectrum in figures 17 
and 18 demonstrates higher sensitivity to the number of spectral groups coupled to 
hydrodynamics. The effective emission radius for the equivalent Planckian flux can be 
evaluated as rem ~ 0.7 mm. Figure 19 shows that it is around this radius that the 
spectral optical depth is on the order of unity. 

Our calculated spectrum in figure 18 appears to be in a fair agreement with the 
observed x-ray spectra for 6 mg cm~^ tungsten arrays tested on the Z machine [11, 4], 
although the published experimental data at hu > 3-4 keV are rather scarce. In fact, 
when we superpose our spectrum in figure 18 on that from [31], we observe a very 
good agreement without even rescaling the absolute fluxes. The experimental points for 
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hu > 2 keV, quoted in [4, 31], do indicate the presence of a hard x-ray component with 
an effective temperature of T^2 ~ 0.6 keV, whereas the main emission is reasonably well 
approximated by a blackbody spectrum with T^i ^ 165 eV [31]. Note that, according 
to our results, particularly in the region hu = 3-6 keV, the spectral slope appears to 
be significantly flattened as compared to the corresponding Planckian fit of the hard 
component — which implies complications for any direct interpretation of the Planckian- 
fit temperature, inferred from the experimental data in this region. 
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Figure 19. (Colour online) Same as figure 16 but for case Z. 
5.4- Calculated ID x-ray images 

Beside spatially integrated emission spectra, of certain interest might be theoretical 
spectral images of the imploding pinch. A selection of such images is shown in figures 20 
and 21 for the time t = 3 ns. Here the radiation intensity Ii, = Iu{s, fi) is plotted as a 
function of distance along an imaginary observation slit, perpendicular to the pinch axis, 
as it would be registered by an observer at infinity; the photon propagation direction 
Q is also perpendicular the pinch axis. Again, these images have been constructed in 
the post-processor mode by separate integration of the transfer equation (7) along a 
predefined set of rays (long characteristics) at selected photon energies. This enabled 
us to get rid of the numerical diffusion inherent in the method of short characteristics. 

Our principal motivation for presenting the images in figures 20 and 21 was to 
illustrate how one could possibly resolve the RD shock front, buried deeply inside the 
imploding plasma column. Figure 20 demonstrates that in the low-mass case A this 
could already be achieved by radiography at photon energies around hu ^ 2 keV. In the 
more massive case Z one has to do the measurements in harder x-rays at hu > 8 keV. 
The softer part of the spectrum reveals only a broad blurred halo from the imploding 
plasma, whose size depends on the observation frequency. 



Radiation- dominated implosion of tungsten plasma 



27 




Figure 20. (Colour online) ID x-ray image of the imploding pinch at two different 
frequencies in case A at i = 3 ns. 
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Figure 21. (Colour online) Same as figure 20 but for case Z and three different 
frequencies. 



6. Summary 

In this work, within the framework of pure radiation hydrodynamics, we have attempted 
to present a detailed physical picture of how the kinetic energy of the imploding high- 
Z (tungsten) plasma in wire array z-pinches could be converted into powerful bursts 
of x-rays. Having concentrated on a self-consistent modelling of the emergent x-ray 
spectra, we adopted the simplest possible formulation of the problem. In particular, we 
assumed that at the final stage of kinetic energy dissipation the dynamic effects due to 
the magnetic field could be neglected, and that the imploding tungsten plasma has a 
perfectly symmetric one-dimensional cylindrical configuration. Both assumptions imply 
severe idealization of the problem, and how realistic are the conclusions reached under 
them, remains to be clarified by future work. 

The reason for our ID statement of the problem is simply because the ID picture 
is always a necessary starting point when exploring a complex physical phenomenon: 
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later on it may serve as a valuable reference case — especially if it manages to capture 
the basic physical features of the studied phenomenon. 

However, even if we skip the initial phase of plasma acceleration by the j x B force 
and stay within the ID picture, there remains a question of possible dynamic and kinetic 
effects due to the (partially) frozen-in magnetic field. We do not expect that such effects 
can significantly alter the present physical picture of the x-ray pulse formation (at least 
not in the phase of what we call the main x-ray pulse) simply because the initial Alfvenic 
Mach number is very high ( > 40). Later on, as the bulk of the imploding mass passes 
through the stagnation shock and the pinch enters the stagnation phase with a Bennet- 
type equilibrium, the effects due to the magnetic field and the ensuing MHD instabilities 
may, of course, become much more significant. This second phase of the x-ray pulse, 
which may in fact account for a large portion of the total emitted x-ray energy and be 
strongly dominated by the MHD effects, was not the topic of our present work. 

Within the approximations made, we have demonstrated that the conversion of 
the implosion energy into quasi-thermal x-rays occurs in a very narrow (sub-micron) 
radiation-dominated shock front, namely, in an RD stagnation shock with a supercritical 
amphtude according to the classification of [9]. We analyzed the structure of the 
stagnation RD shock by using two independent radiation-hydrodynamics codes, and 
by constructing an approximate analytical model. 

We have found that the x-ray spectrum, calculated with the 2D RALEF code by 
solving the equation of spectral radiative transfer in the imploding plasma, agrees fairly 
well with the published experimental data for the 6 mg cm~^ tungsten wire arrays tested 
at Sandia. The hard component of the x-ray spectrum with a blackbody temperature of 
Ti.2 ~ 0.5-0.6 keV is shown to originate from a narrow peak of the electron temperature 
inside the RD stagnation shock. Our approximate model clarifies how the width and 
the amplitude of this temperature peak depend on the imploding plasma parameters. 
The softer main component of the x-ray pulse is generated in an extended halo around 
the stagnation shock, where the primary emission from the shock front is absorbed and 
reemitted by the outer layers of the imploding plasma. 

In reality, due to ffow non-uniformities, the narrow front of the stagnation shock 
will almost certainly have a much more irregular shape than in the present ID picture. 
But its main characteristics — the transverse thickness Avg and the peak electron 
temperature — are controlled by the plasma ffow parameters (the implosion velocity 
Uo, the mass ffux density PqUq, the plasma thermal conductivity k and the spectral 
absorption coefficient k,^) that are not expected to be dramatically affected by moderate 
ffow perturbations. Hence, we expect that radiation-hydrodynamics simulations of 
realistically perturbed implosions should produce emergent x-ray spectra close to those 
calculated in the present work. 
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